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Abstract 


We outline a model of abelian-Higgs strings with variable scalar and vector core radii. 
In general, the functions determining the time and position-dependent core widths may 
be expressed as arbitrary left or right movers, of which the usual constant values are a 
particular solution. In this case the string may carry momentum, even if the embedding of 
its central axis remains fixed, and the resulting objects resemble “necklaces”. Some possible 
astrophysical applications of lumpy strings, including as potential engines for anomalous 
gamma ray bursts, are also discussed. 


1 Introduction 


In 1973, Nielsen and Olesen’s seminal paper, demonstrating the existence of a topologically 
stabilized vortex solution in the abelian-Higgs model [T], marked the beginning of the theory of 
cosmic strings as topological defects, formed during symmetry breaking phase transitions in the 
early universe mm- Since then, an enormous number of string-type solutions have been found 
in many different field-theoretic models, including electroweak strings in the Weinberg-Salam 
model H El E], and strings in GUT scale theories [7]. Phenomenologically, superconducting 
strings are also of great interest in cosmology and astrophysics BM- 

Interest in cosmic strings was also renewed after the discovery that T-strings (fundamental 
strings governed by the Nambu-Goto action) and D-strings (one-dimensional D-branes), could 
be formed in large numbers at the end of inflation, in string theory inspired models of the early 
universe Though not field-theoretic, such strings are nonetheless “cosmic” in origin and 

provide an alternative paradigm in which to study the cosmological implications of the existence 
of string-like objects. 

With the exception of “necklaces” m-m which, in many models, are a form of hybrid 
defect; that is, a series of monopole-like “beads” connected by strings, rather than strings in 
their own right, virtually all these models have one thing in common. Whatever the underlying 
nature of the string, it is, essentially, a “one-dimensional” object. Though, strictly, this is only 
true for cosmic superstrings, at the classical level, while topological defect strings necessarily 
occupy a three-dimensional volume, in the latter the string radius is assumed to be constant and 
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much smaller than the string length, so that it may be neglected on cosmological scales. Even 
in the case of necklaces, the monopoles do not typically occupy large volumes, or have widths 
much larger than the string itself. 

An alternative model was proposed in [20], in which a defect string with varying local 
radii was found to exist. By modifying the standard abelian-Higgs model to include spatially- 
dependent couplings for the scalar and vector fields, static, non-cylindrically symmetric solutions 
of the resulting field equations were obtained. Though phenomenologically interesting, in the 
absence of any conclusive proof for the spatial (or temporal) variation of fundamental constants 
[211 [221 [23], such solutions must be considered as toy models. 

In this paper, we propose a new solution of the standard abelian-Higgs field equations, which 
represents a generalization of the Nielsen-Olesen string. For a string whose central axis lies along 
the coordinate z-axis, this solution is characterized by the addition of two arbitrary functions 
of z ± f, which represent arbitrary left or right movers, and which determine the effective local 
radii of the scalar and vector cores. The generalized solution thus represents a “lumpy” cosmic 
string, and it is found that the lumps propagate at the speed of light. Phenomenologically, the 
solution combines several features of existing models in an interesting way. Specifically, it is 
non-cylindrically symmetric, like the static solution presented in [20|, but possesses momentum, 
as well as mass-energy, which may be interpreted as a kind of neutral current, as in supercon¬ 
ducting strings. In addition, for lumps that are highly localized along the string length, it may 
resemble a necklace, with highly energetic traveling beads. 

Importantly, however, unlike necklace models, there is no theoretical maximum value for the 
radius of the lumps, at least not in the classical theory considered here. In principle, quantum 
effects may limit its value to the same order of magnitude as the standard string width, though 
further detailed analysis is required to establish whether this indeed the case. Therefore, in order 
for lumpy strings with conserved momentum to form, all that is required is for local variations in 
the vortex radii, at the epoch of string formation, to be correlated, rather than random, at least 
over some characteristic length scale. Again, only a detailed and, possibly, quantum mechani¬ 
cal analysis of the field evolution during the symmetry breaking phase transition can establish 
whether this is possible on scales that will effect the macroscopic string dynamics. However, 
we assume, for now, that such correlations can exist, and determine the general features of the 
resulting lumpy string solution. 

The structure of this paper is then as follows. In Secs. 12.11 and 12.21 we review the abelian- 
Higgs model and the Nielsen-Olesen solution, respectively, while in Sec. 13.11 we present the 
lumpy string ansatz and demonstrate that the field equations may be solved consistently. In 
Sec. 13.21 we determine the lumpy string energy and momentum, as well as the integrated values 
of the local pressures and shears. This allows us to estimate the energy and momentum emis¬ 
sion as the string relaxes to the (lower energy) Nielsen-Olesen configuration, in Sec. 14.H and 
some possible astrophysical implications of this are considered in Sec. 14.21 In particular, we 
consider the possibilities that energy emission from lumpy strings may provide an explanation 
for anomalous gamma ray bursts [24] . or contribute to the observed positron excess coming from 
the galactic centre [25] Finally, a brief discussion and summary of our results is given in Sec. [5j 
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2 Background: The abelian-Higgs model and the Nielsen-Olesen 
string 

2.1 Recap of the abelian-Higgs model 


In natural units, and using the metric convention (H-), the abelian-Higgs action is 

S = I d^x^g - VM)^ , (2.1) 

where G {0,1,2,3} and H(|(^|) is the potential term, which is determined by the symmetry 
breaking energy scale, rj, and the scalar coupling, A: 

= ( 2 . 2 ) 

Following the conventions of Vilenkin and Shellard [2], the gauge covariant derivative and 
the electromagnetic field tensor F^y, are defined according to 

D^ = d^- ieAf,, F^y = d^Ay - dyA^, (2.3) 

where e is the vector coupling, and the covariant equations of motion (EOM) are 

+ Ac/) (|(/)|2 - r?2) = 0, (2.4) 


du {V^F>^'^) + m - 4>D^(tA) = 0. 


The action, Eq. ()2.1I) . is invariant under local U{1) transformations of the form 




A^ —)• A'^ — A^ + —dfj,A{x), 


where A(x) is any single-valued real function, giving rise to the conserved current 

= -ie - (l)D^J) . 


(2.5) 

( 2 . 6 ) 

(2.7) 


The conservation of is expressed directly via the vector EOM, Eq. (USD, and the corresponding 
conserved charge is 




( 2 . 8 ) 


The energy-momentum tensor is defined implicitly by varying the action with respect to the 
metric, 


6S = j T>^^5g^,y/^d^x, (2.9) 

so that 

^ d{Cy^) ^ - g^^''C, ( 2 . 10 ) 

c>9t^v 

where C is the Lagrangian density; that is, the sum of terms inside the curly brackets in Eq. 

(EU). 
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2.2 Recap of the Nielsen-Olesen solution 


In cylindrical polar coordinates (t, r, 9, z) and assuming a Minkowski background so that, 

ds^ = rj^ydx^dx'^ = c^dt^ — dr"^ — r^d9^ — dz^, ( 2 - 11 ) 

the ansatz for the Nielsen-Olesen string is [T] 

(/)(r, 9) = r?/(r)e*"®, Aq = -ae{r), Ar = = Aq = 0, (2.12) 

e 

where /(r) and ag(r) are dimensionless functions which obey the boundary conditions 

0, r = 0 


= 


ag{r) = 


1, r 
0, r = 0 


oo, 


1, r —>■ oo. 

These give rise to the specihc form of the EOM [2] 


dR^i I dr 2 \ r. 


''s|n| 

d?ag 

dW~, 

t;|n| 


1 


s|n| 


dag _j_ 2 ^ '^v\n\ 




^v\n\ dRy\n\ 


f\l-ag) = 0, 


where we have defined 


^i\n\ 


^i\n\ 


i e {s,u} , 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 


and where r^i^j and r.u|„| denote the values of the scalar and vector core radii for an |n|-vortex 
string, respectively. We also use the simplified notation, and ry, to refer to the scalar and 
vector core radii of |n| = 1 strings and define the parameters 


P\n\ 


2 

Un\ 

o 1 P 9 ; 

r I I ri 

s|n| s 


(2.18) 


for later convenience. Eqs. (j2.15l) - (j2.16p are solved, to leading order in the uncoupled regime, 

by [51 [TO] 


/(r) « ■ 

1 {r/rs\r. 

I'jhI 

i\} 5 

1, 

r < rs\n\ 

^ ^ '^s\n \) 

(2.19) 

ag{r) « 

1 {r/ry 

h)^ 

1, 

r < Vyi^i 

r > r^|n|- 

(2.20) 


Strictly speaking, the power law solutions for /(r) and ag{r) given above are only valid in the 
ranges r <C i G {SjU}, while the asymptotic forms hold only for r S> rj|„|. However, 

by assuming that each holds, approximately, up to the limiting value r = rj|„|, and ignoring 
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the discontinuity of the derivatives, we can easily obtain order of magnitude estimates for the 
physical parameters of the string without the need for detailed numerical calculations. 

The topological winding number, n G Z, is obtained as 

1 d'& 

"=2^X 

where i? represents the phase of the scalar field (j). In general, the core radii are expected to 
be of the order of the Compton wavelengths of the associated bosons, and m~^, but may 
contain some |n|-dependence. This dependence is expected to take a simple form [I9l[26l[27] so 
that here we assume 


rs\n\ ~ |n|^m/= |n|^(VAr/) \ 

rv\n\ ~ \n\"rn~^ = |n|"(ery)"\ (2.22) 

where ^ > 0, e > 0 are constants. 

Plugging the ansatz Eq. (I2.12p . together with the approximate solutions in Eqs. (12.19^ - 
(j2.20l) . into Eq. (I2.10p . the only nonzero component of the energy-momentum tensor is 
so that the only constant of motion is the Hamiltonian, E|„|, which depends on the absolute 
value of the topological winding number. The (constant) mass-energy per unit length of the 
string, = f T^^rdrdO, is simply the total energy divided by the total length and, for 
and ry^n\ given by Eq. (I2.22p . it is straight forward to show that, assuming (/^ > 1); 

which corresponds to a Type II superconducting regime [2], 


H\n\ « 27rr/^ 


n -|- 


I |2-2e ,1 |2f , 

n -I- n ^ -|- 


nP In 


P\n\ 


(2.23) 


However, in pioneering work [261 [27], Bogomol’nyi showed that, for a Nielsen-Olesen string at 
critical coupling (e = y/X, /3 = 1), //|„| must satisfy the condition 

M|n| > (2.24) 

otherwise known as the Bogomol’nyi-Prasad-Sommerfield, or BPS bound [28]. Since stability 
implies saturation of the bound [2] (see also [29] for further work on the stability of abelian-Higgs 
strings), we require 


0<?<l/2, e>l/2, (2.25) 

for \n\ ~ 0(1). Alternatively, we may set 

C = e = 1/2, (2.26) 

for arbitrary |n|. It is common in the literature to take <7 = 0 and e = 1/2, following Bogo- 
mol’nyi’s original assumptions [261 [27], which satisfy the conditions in Eq. ()2.25p . but in the 
present work we adopt the latter conditions, given in Eq. (|2.26p . to allow for large |n|. For any 
applicable choice of parameters, the order of magnitude estimate for the mass-energy per unit 
length of a Nielsen-Olesen string with winding number n, at critical coupling, is 

M|n| ~ 27r7/^|n|. (2.27) 
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Though it is not physically necessary for the components of the total tension vector within the 
string = f T^j.fdrdO, T® = f T^grdrdO and = f T^z'f'drdO (also in units of energy per 
unit length), to be conserved, is straightforward to show that, for the Nielsen-Olesen string 

fi\n\ = -T", (2.28) 


for any value of /3|„|, and 


V = T^ 


— 27T'lf 



2 - |npln (yMh) 


(2.29) 


Hence, for ^ = e = 1/2, as in Eq. (|2.26p . and at critical coupling, the only nonzero component 
of the tension is which is equal to minus the mass-energy per unit length of the string. 
Substituting the ansatz Eq. (I2.12|) into Eq. (12.7p . it is also straightforward to show that 
= 0, V// G {0,1,2, 3} and hence that the string is uncharged. 


3 Field configuration and constants of motion for a lumpy string 


3.1 Field configuration for a lumpy string 


In this section, we consider a model of an abelian-Higgs string, with variable scalar and vector 
core radii, in which the essential phase structure of the Nielsen-Olesen solution is preserved. 
That is, we wish to preserve the circular symmetry and topological winding number of the 
string cross section at all times and at each point along the string length, even if the radius of 
the cross section is itself a function of t and z. We therefore take the ansatz 

^{r,e,z,t) = 

(3.1) 

Tl 

Aq = -ag{r, z, t), Ar = Az = Aq = 0. 
e 

(3.2) 

The scalar EOM then becomes 

IdF n'^F^ ^ 

dr2 ^ r dr r2 ^ + dz'^ dF 2r1 ^ ^ ~ ' 

(3.3) 

and the 0-component of the vector EOM is 

d'^ae Idag d'^ag d‘^ag 2^2 

Q,p2 ^ Qj, Q^2 Q-f-2 j.2 

(3.4) 

The z- and t-components both give 

d'^ag d^ag _ 
dz'^ dF ’ 

(3.5) 


which allows us to cancel these terms from the 0-component, and the r-component is identically 
satisfied. 
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Clearly, there exist functional forms for F{r, z, t) and ag(r, z, t), which satisfy the EOM, that 
include both left and right movers, i.e. 

F{r,z,t) = F{r - gs{z±t)), ae{r, z,t) = agir - gy{z ±t)), (3.6) 

where gs{z ± t) and g^{z zh t) are arbitrary functions oi z F t and which, in principle, may 
include superpositions of the form gs{z F t) = gs{L){z F t) F gs(R){^ ~ t)- However, we must 
carefully examine the physical nature of these waves. In Cartesian coordinates {t,x,y,z), an 
ansatz representing waves propagating on the string, that is, representing displacements of the 
string central axis from the original coordinate z-axis, takes the form [30] 

(p{x,y,z,t) = Af,{x,y,z,t) = Ai,{X,Y), (3.7) 

where 


X = x - gxizFt), Y = y - gy{zFt), (3.8) 

and <I>(x,y), Ag{x,y) represent the original Nielsen-Olesen solution, expressed in Cartesians. 
Such a solution represents a string, of finite width, whose central axis is located at the “origin” 
with respect to the transformed coordinates X = Y = 0, rather than at the origin of the physical 
coordinates of the background space, x = y = 0. In this case, the general string configuration, 
at a given time t, is that of a curved cylinder, but each vortex cross section (taken perpendicular 
to the central axis) remains a circle of fixed radius. Solutions of this type, known as “traveling 
wave” solutions, for abelian-Higgs strings, were shown to exist in [3U] (see also the appendix in 
[31] for further discussion). In fact, to obtain true traveling wave solutions, it is not sufficient 
to have only nonzero (j), Ax and Ay, but nonzero A^ and of the form 

^0 = Mdy^xiX, Y) F gxAy{X, H)] = FA^, (3.9) 


are also required. In this case, though the functions gx{zFt) and gy{zFt) may represent either 
left or right movers, the choice of sign in the argument of each must match that in Eq. (13.9j) . 
so that linear superpositions of waves propagating in opposite directions are not allowed m- 
As we shall see, similar restrictions apply to the left/right movers that characterize the z- 
and t-dependence of the lumpy string ansatz but, first, let us take a closer look at the solutions 
proposed in Eq. (13.61) . Since, by analogy with Eqs. (|3.7p - (l3.8l) . the values of gs{zFt) and g^izFt) 
represent the displacement from r of the variables TZg = r — gs{z F t) and TZ^ = r — gx{z F t), 
respectively, at the point z at time t, we see that the conditions TZg = 77„ = 0 do not correspond 
to a unique point within a cross section of the string. Moreover, since the true vacuum can exist 
within a vortex at only one point [2], we are required to impose the boundary conditions 


F{r,z,t) = I J 


r = 0 

r —>■ oo, Vz, t 


(3.10) 


ae{r,z,t) 


0, r = 0 

1, r —^ oo, Vz, t, 


(3.11) 


which is impossible, for the solutions proposed in Eq. (13.6p . without setting gs = Qv = 0, Vz,t. 
Such solutions, though they exist mathematically, therefore violate the physical requirement 
that {(p) = 0 at one point only within a given vortex [2]. 

As an alternative to Eq. (13161) . we may instead search for approximate solutions of the form 


F{r, z,t)= F t)=aei R' 




(3.12) 
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where we have defined 



r 



i G {s,z;} , 


(3.13) 


and where and represent the effective local values of the scalar and vector 

core radii (for given values of 2 : and t), respectively. The physical nature of the waves described 
by the ansatz, Eqs. (ISID-dMl), now becomes clear and we may write r^j^j( 2 ;,t) and in 

the form 


= '^i\n\ 9 i{z, t), ie {s, u} , 


i\n\ 


(3.14) 


where the gi(z,t) are again arbitrary (though this time, dimensionless) functions. Substituting 
for F{r,z,t) and a 0 {r,z,t) from Eq. (13.121) into Eqs. (I3.3I) - (I3.5I) . the scalar EOM becomes 


d^F 




+ 


1 dF 


rrF 


il-aeY- 


I ^ff 
1 / ^ I I 

2 I s n 


RllIdRil {Rl\n\)^ - -- 2 


F(F2 - 1) 




e //\2 


+ (Kill) 


d?F 




'K\IX (K\n\\ 

dz I 1 I 


- R 


e// 


dF 


eff 


5 |n| s\n\ 






- 2 


dr'^f { 

dz 


dr^f { \ 

dt I 


2-1 




= 0 , 


and the 2 - and t- components of the vector EOM are equivalent to 

2 / ,rr\ 2n 


/pe//x2 ^ 




- R 


eff 


dag 


eff 


dr^dd 

v\n\ 

dz 

ref 
v\n\ 


dr^f d. \ 

■eH I 

dt I 


q 2 eff q 2 eff ' 
..i„i v\n\ 


dR^df I 1 dz'^ 

^v\n\ ' ' 


dF 




- 2 


dr 


eff' 

i»|n| 


dr 


eff' 

i;|rx| 


dz 


dt 


= 0 , 


while the 0-component takes a form analogous to Eq. (I2.16p . under the identifications ag ag, 
f ^ F, rllJ^^{z,t) ^ r^\n\ and t) o Ry\n\- 

The EOM therefore reduce to 


d^F 


d{R^df ^2 


+ 


1 dF 


tFF 


*'s|n| ^ 


RildRil {R^df \2 


s|n| 




{l-aoY-- 


, eff ' 
1 f I \ 

2 I 


F{F^-l)=0, (3.15) 


d?ag 


dag 


rYI dR^d, 

v\n\ v\n\ 


+ 2 



2 


FYl-ag)=0, 


if the functions gi{z,t) satisfy the conditions 


% _ I % 

\dz) \dt J ’ dz dt' 
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(3.16) 


(3.17) 










































which, in turn, automatically imply 


^ d'^gj 
dz^ dt'^ 


(3.18) 


However, Eq. (13.171) is the stronger condition, since Eq. (I3.18P is satisfied by superpositions of 
arbitrary left and right movers (i.e. functions of the form gi{z,t) = gi[L){z + t) + gi[R){z — t)), 
whereas Eq. (|3.18l) is satisfied only by arbitrary left or right movers: gi{z,t) = gi{z + t) or 
gi{z,t) = gi{z -t). 

Clearly, by analogy with the solutions obtained for Nielsen-Olesen strings, Eqs. (I3.15p - (I3.16I) 
are solved, to leading order in the uncoupled regime, by 




F{r,z,t) 


> r^f 




(3.19) 


1. r ><"(.,*). 


(3.20) 


Physically, these functions represent the spatial and temporal variation of the scalar and vector 
core radii, for a string with a fixed central axis (relative to the background space), lying parallel 
to the z-axis. Combining the conditions in Eqs. p3.17p - p3.18l) with the solutions given in Eqs. 
(I3.19p - (I3.20I) automatically implies 


dF 

dz 


'm) ’ 


dF _ dF 


(3.21) 


dae 

dz 


and, though it follows that 


QaeV dae _ ^ dap 


dt 1 ' dz 


d2F _ d^F 
dz"^ dF ’ 

di^ae d'^ag 
dz"^ dF ’ 


dt ’ 


(3.22) 

(3.23) 

(3.24) 


Eqs. (j3.21l) - (j3.22l) represent stronger constraints, which prevent the possibility of superpositions 
of left and right-movers in the radial oscillations, as stated above. 

Einally, before moving on to calculate the energy and conserved momentum of the lumpy 
string in the next section, we note two further restrictions that physical considerations place on 
the mathematical form of the functions gfizfi). Firstly, since we have used an ansatz describing 
(formally) infinite strings, we must impose periodic boundary conditions in the z-coordinate 
and consider a finite section of length A when determining the constants of motion, in order to 
prevent the energy and momentum from diverging for 2 ^ ±oo. Similarly, the fact that both 
quantities must be hnite and conserved, for string sections of finite length, implies periodicity 
in t in order to prevent divergence for t —)• oo. We therefore impose the conditions 


gi{{z -4 mXfij ± {t+pT)) = gfiz ± t), m,p G Z 


(3.25) 
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for some arbitrary length scale and time period T, which may be associated with a frequency 
via 

27r , 27r , , 

uJz = — = ±^2 = ±-p. (3.26) 

1 Az 

Secondly, the finiteness of the energy density also requires us to choose the form of gi{z±t) such 
that 

giiz±t)^0, 'iz,t (3.27) 

as we shall see explicitly in Sect. 13.21 In the classical theory therefore, an arbitrary cut-off must 
be imposed by hand but, thanks to the well known quantum properties of abelian-Higgs strings, 
we may simply set 

gi{z±t)>l, yz,t (3.28) 

which ensures that r^^^J^{z,t) > '^z,t. This ensures that the local radii of the scalar and 

vector string cores do not drop below the minimum values dictated by the Compton wavelengths 
of the associated bosons. For the specific choice of gs{z ±t) = gy{z ± t) = 1, we recover the 
Nielsen-Olesen solution [T]. 


3.2 Lumpy string energy and momentnm 


If a = {0,1, 2, ...d}, denotes a set of D = d + I Killing vectors, then 

denotes a set of D conserved currents, 

= 0 , 


(3.29) 


(3.30) 


and the associated conserved charges are 

n“ = y = J (3.31) 

For the sake of completeness and clarity, we list all independent components of the energy- 
momentum tensor below: 







(1 — og)^ -|- 


4^2 





If 


,0r ^ _ 2 I 

^ \ dr dt dr dt ] ' 

TOd ^ 


(3.32) 

(3.33) 

(3.34) 
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= - 77 ^ <j 2 


2 j ^dF dF |npr^ daQ dag 

dz dt dz dt 


(3.35) 


^ U 5r 




(1 — oe)^ + 


3 Inpr^ /(9a6i\ ^ 1 

4 ^2 \ dr J Ar'i 


(F^ - 1)^ 


(3.36) 


T'’® = 0, 


(3.37) 


2 I ^9F 5F |npr^ dag dag 


^ (9r 02 ; dr dz 1 ^ 


(3.38) 


= 7?^ < - 


1 fdFy |npF2 


^2 \ 


H-- 4—(1 - ae) + 


3 InPr,^ /5oe 


- 


4 ^4 \ dr J 4r^r^ 


r(F2-l)2^, (3.39) 


F^^ = 0, 


(3.40) 


^ (f )'^ ^ (1?)^ - if -') 


(3.41) 

For (3 + l)-dimensional Minkowski space there are ten Killing vectors, namely four related to 
space-time translational invariance, three to spatial rotational invariance, and three to boost 
invariance, respectively. However, from Eqs. (|3.32p - (l3.4ip . it can be seen that the only relevant 
independent Killing vectors, giving rise to nonzero conserved quantities for the “straight” lumpy 
string described by Eqs. (I3l|)-(I321) , are 

= (1,0,0,0), = (0,0,0,1). (3.42) 


Thus, the independent conserved quantities are the Hamiltonian F|„| and the momentum in 
the ^-direction, = —Pz- The total 3-momentum is F|„| = y—P^Pi = y/—P^Pz, and the 
conserved 4-momentum is given by 



We may also define the integrated tensions (negative pressures) as 


(3.43) 


T^ = j T\i^y^d^x, 


(3.44) 


where the notation T*(q indicates that no sum is implied by the repeated index, and the inte¬ 
grated shears between two arbitrary directions, x* and , via 


Sd 


= j T^y^d^x, 


(3.45) 
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for i,j G {r, 6 , z}. 

Using the approximate solutions for F{r,z,t) and ag{r,z,t), Eqs. (l3.19p - (|3.2Up . the order of 
magnitude values for the nontrivial constants of motion, integrated pressures and shears for a 
finite section of lumpy string, of length A, are 




S 


rz 



dg^ \ 
rv\n\9l dz J 


= 0 , 


(3.50) 


where we have used the approximation |n| +1 ~ |re| for |n| > 1, and ignored numerical factors of 
order unity multiplying individual terms within the integrands. The vanishing of the integrated 
shear follows directly from the imposition of periodic boundary conditions on gs{z,t) and 
gy{z,t), Eq. (I3.25p . if we choose to set 


A = mXz, (3.51) 

though a nonzero local shear, S^^{z,t) = f T^^{z,t)rdr, is, in general, present. 

It is straightforward to check that the standard expressions for the constants of motion and 
integrated pressures/shears for the Nielsen-Olesen string are recovered for gs{z,t) = gc{z,t) = 1. 
However, Eqs. (|3.46l) - (|3.50p also simplify considerably under the more general assumption 


9 s{z,t) = gv{z,t) = gc{z,t), 


(3.52) 


together with critical coupling and ^ = e. In this case. 


eff/ 

"s|nV 


(/^i = 1) 


eff 
— f JJ 
c\n\ 


|n| 


(3.53) 


and the logarithmic terms vanish, while those involving gg, gv, ?’s|n| a'lid can be consolidated. 
For the specific choice, ^ = e = 1/2 (as before), and using the fact that dgddt = Fdgddz, we 
then have 


E 


\n\ 


27rr/^|n| A 



\n\9l + ^ + 


c\n\ 



(3.54) 
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(3.55) 


« ±27rT/^|n|r^l„l / dz 
Jo 




= T 


e 




(3.56) 


« 27rr?^|n|A 

where the choice of sign in Eq. (|3.55l) is positive for right-movers and negative for left-movers. 
Given also that gc{z,t) > 1, we may, to lowest order in the approximation, neglect terms in g~^ 
and g~‘^{dgc/dz)‘^, as long as their signs do not conflict with those of terms proportional to g^ 
and g'l{dgcldzY, respectively. 

Furthermore, by rewriting the 2 ; coordinate in terms of a dimensionless parameter a € [0, 27r), 
such that 


1 + 


f 


dz 


1,2 1 

Hdc + — 


- 

^c\n\ 


dgc\ 
dz ) 


1 + 


(3.57) 


z(a) = (27r) ^Afj, 


(3.58) 


we may write 


gdz ±t) = gdkzZ -|- uj^t) = gd'ma + Uzt), 


(3.59) 


where ooz = J^kz = ±27r/A2 = ±27rm/A, in accordance with Eqs. (I3.26P and (I3.5ip . Of course, 
different plane wave modes in the Fourier expansion of gdzH) correspond to different frequencies 
ujj and associated wavelengths Xj, such that Uj = ±kj = ±2TTjXj = :iz 2 Trmj jA, for some 
rrij G N. By writing gdz ± t) in the form given in Eq. (|3.59D . which implies dgc/dt = ojzdgc/du, 
dgcjdz = kzdgc/du, dgdda = mdgc/du, where u = kzZ + (jJzt = ma + cozt, we select a single 
mode as being “characteristic” of the wave form. The most natural choice is the mode with 
the highest amplitude, whose wavelength will give the approximate width of any “lump” on the 
string. The integer m may then be interpreted as defining the number of lumps present in a 
string section of length A. This gives 


El 


|n| 


2'Krf\n\A 


1 + 


1 

27r Jo 


Hal + 


(2vr)^ 


c\n\ 


A2 


{daacf \ da 


(3.60) 


K. ±27r77^|n|A x 




A 2 2-k 


^ 5“ / {d^gd) da, 


(3.61) 


= T^ K -2TTgdn\A x 


iCPi 


da, 


(3.62) 


Ri — 27r7/^|n|A 
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r2TT 


(27r)2r2 

A^ 




-(^affc)^ > da 


(3.63) 
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We may then write 

2 p27t 2 ^27r 2 p27r 

^ Jo ^ ^ i) ^ ^ Jo ^ 

where o? > 1 , 7 ^ > 0 , <5^ < 1 are constants, with equality holding if and only if gc{z,t) = 1 , 
Vz, t, so that 


E\n\ ~ 27r7^|n|A 


, , o (2vr)V5„|72' 
1 + \n\a^ + ■ 


A? 


« ±27r7^|n|A x 


(2-)^r5„|7^ 

Ai 


= T^ ^ -27r7^|n|A(a^ - 5'^), 


K, —27r7^|n|A 


1 + |n|a"^ — 


2 (2-)^in|7" 


A 2 


(3.65) 


(3.66) 

(3.67) 

(3.68) 


and 


n 


|n| 


27r7^|n|A(l + |n|a^) 



(27r: 


i2.^2 


' I I 

c\n\ 


A? 


27^ Y 
1 + \n\a‘^ J ' 


(3.69) 


Interestingly, it is possible for the integrated pressure in the 2 ;-direction to be zero, even in the 
nonlinear case, if 


A^ = ( 2 vr) 


2 2 

^c|n| 


7 


1 + n W 


(3.70) 


but can only vanish for the Nielsen-Olesen string. However, in general, there are 

local/instantaneous momentum densities, tensions and shears, which depend on both a oc z and 
t, i.e. 


g\n\{u) « 27rr]‘^\n 


l + \n\gl{u) + 


( 2 ^)' 


c\n\ 


A^ 



(3.71) 


,0 2 , ^c\n\ (dg, 

V {u)^±2.g\n\x^^[^—j , 


(u) = T^{u) ^ —27r?7^|n| 


ffcH - 


9l{u)\ 


(3.72) 

(3.73) 


T^{u) « —27r7^|n| 


( 2 ^) ^c|n| ^ 

du 


1 , I I 2/ \ ' ' c|n| I u.yc 

1 + \n\g^{u) - 


AI 


(3.74) 
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(3.75) 


where 




-(27r)Wkc|n| 


dgc 
du ’ 


u{t,a) = ma + ojzt. (3.76) 

In general, therefore, we find that 

g\n\{u) > \r^{u)l (3.77) 

where equality holds only in the Nielsen-Olesen case, for which both the tension and the mass- 
energy per unit length become constant. From the /i = 0 component of the vector EOM, it is 
also straightforward to see that, for the lumpy string, = 0 , and the string is uncharged. 


4 Cosmological and astrophysical implications of lumpy strings 


In this section, we consider some possible astrophysical and cosmological implications of lumpy 
cosmic strings. As a first step, we give a simple estimate of the energy and the momentum 
emission for the lumpy string conhguration, as it relaxes to the standard Nielsen-Olesen form. 


4.1 Energy and momentum loss of lumpy cosmic strings 


Under the simplifying assumptions, gs = gv = Oc, ^ , C = ^ = 1/2, made in the previous 

section, it is straightforward to estimate the loss in energy AEj^i and momentum AP|„|, over a 
string section of length A, as the string relaxes to a Nielsen-Olesen configuration. In terms of 
the new model parameters a^, 7 ^ and A^, these are 


AE 


|n| 


27rr7^|npA 


+ 


(2vr)V,V ^ 

A2 



(4.1) 


APi 


\n\ 


2'Krf‘\n\^A (27r)^rc7^ 

^ a! 



(4.2) 


where g is dimensionless, {G/c^)~^ = mpc^/lp and, from now on, we write fundamental con¬ 
stants such as G, c and h explicitly. An estimate for the characteristic time scale for the loss 
may be obtained from the characteristic frequency, 


At « T 


27r 

LOz 


27r 

kz 



(4.3) 


Roughly speaking, characterizes the additional rest mass-energy of the lumpy string 
(over and above that usually contained in a length A of Nielsen-Olesen string), while 7 ^rg/A^ 
characterizes the additional kinetic energy due to the propagation of the lumps. The latter 
is higher when the average distance between lumps, Az, is small, and the average square of 
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the gradient, 7 ^ is large. In other words, when there are many, closely spaced, sharp peaks 
in the radial profile. Alternatively, when the lumps are widely separated but highly localized, 
corresponding to a sharp peak in over a flat background, it would be more accurate to 
introduce a two-scale model, in which the characteristic inter-lump distance and the average 
lump width appear. 

Nonetheless, in our current one-scale model, the difference between these two scenarios simply 
corresponds to a smaller or a larger value of 7 ^, respectively, for a given value of A^. Therefore, 
for any model in which the lumps are characterized by sharp peaks, it seems reasonable to 
assume that the ratio is large, even if we take > Tc as a fundamental cut-off due to 

quantum effects, and that the kinetic term is of leading order. For radial profiles that are slowly 
varying, either 7 ^ is small and A^ is large (relative to Vc), or A^ is small and 7 ^ is very small, 
corresponding to extremely shallow ripples on the string. In either case, the ratio 7 ^rg/A^ 
remains small, so that the rest mass term o? dominates. We therefore investigate these two 
limiting regimes separately. Hence, 


AH 


\n\ 


27r?7^|n|^A x 


(2vr)Vg7^ 

A2 


X 



< 


(2vr)Vh^ 

A2 


(4.4) 


AH|„| 27r?72|np2Aa2 X ^ 

In the first case, we have 

AH|„| « AP|„|C, (4.6) 

so the emitted energy is mostly in the form of “radiation”, i.e. either gravitational or electromag¬ 
netic waves, or massive but highly relativistic particles. At present, it is unclear whether gravita¬ 
tional radiation (see mm and references therein) or gauge particle emission, which then sources 
metric perturbations [321I33| . is the primary energy loss mechanism for oscillating Nielsen-Olesen 
strings. It probable that, whichever mechanism dominates for traveling waves or oscillating loops 
of cylindrical string, the same mechanism dominates for lumpy strings in this regime. However, 
the opposite limit, Eq. (|4.5p . we may assume that most of the energy is channeled into the 
production of (slow moving) massive particles, so that 

AHj^i Ri AA/|„|C^, (4.7) 

where AM|„| is the total mass of the particles produced, which also depends on the absolute 
value of the winding number, |n|. 

Explicitly writing the correct fundamental constants, the expression for in terms of 
and is 

~ (^) (^oeN^)"^ = , (4.8) 

where the Planck length Ip is defined as Ip = hG/c^, and = e^/qp, where the Planck charge 
is Qp = dvreohc = Aph/cgL^. If e denotes the standard charge of the electron, the dimensionless 
quantitity e,. is equal to the fine structure constant, ag ~ 1/137, though, in principle, we can 
consider toy models in which the coupling is arbitrary. 

We may then estimate the energy emitted by the string section of length A, in either scenario 
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above, in terms of and the appropriate numerical constant, 7 ^ or o?^ respectively. Using 
A = Ax, this gives 


(A£'|„|),.£ici ~ ‘I'K 


. he 


I i2 2 

|n| 7 m 


{^E\n\) 


1 he 


I i2 2 

n a m, 


matter 


2 el 


The corresponding power outputs are simply 


(4.9) 


{^E\n\)rad ^ (AU|„|)matter 1 |npa^m t a 

(At) XI ’ (At) ^ ^ 

In matter dominated emission, the total energy release is proportional to the lump width 
whereas, for radiation, it is inversely proportional, as expected intuitively. As a result, the 
expected power output (per lump) for the former tends to a constant value, whereas, for the 
latter, it is inversely proportional to the lump width squared. 


4.2 Astrophysical applications 


Considering a single lump (setting m = 1 ) and assuming that e = 1.602 x 10“® C, the standard 
charge of the electron, the radiation and matter-dominated power outputs from a string of length 
Ax can be written in an equivalent form as 


')rad 

Af 


8.906 X 10®^ X |n7^ x 



(4.11) 


( AUj^l )mQ,j4ej. 

Af 


2.255 X 10®^ X |n|7^ x 



(4.12) 


It is interesting to compare the above estimates with the energy emission of superconducting 
cosmic strings. Superconducting strings, moving through a cosmic plasma in the presence of a 
magnetic field B, are sources of synchrotron radiation, which may be observed at centimeter 
wavelengths. The energy of this radiation can be estimated using the formula [HI [341 (35] 


AUss 

At 


= 2 X 10 ^^ erg/s. 


(4.13) 


where R 20 = R x 10^® cm is the radius of the string loop, i ?_6 = Bj (l0“® G), and Up is the 
number density of the cosmic plasma particles. For a typical vacuum situation, B-q = 10®, 
R 20 = 1 and Up = 10® cm“®, giving AEss = 2 x 10^® erg/s. For the case of a very dilute 
medium, with R 20 = 3, i ?_6 = 1 and Up = 10 cm“®, we obtain AEss ~ 10®^ erg/s. The 
radio emission of superconducting strings has a very peculiar structure, in that they should be 
observed as a closed line. 

Thus, for |np ~ 7 ^ ~ 0(1), the emitted power of a single-lump string of length Ax = 10^^/p, 
decaying via radiative emission, is greater than or equal to that of a superconducting string 
of equal length. This implies Ax ~ 10^® cm, and the energy is emitted on a timescale At < 1 
millisecond, after the epoch of string formation. However, clearly, our original assumption that 
radiative emission dominates is invalid for Ax Vc, and the latter will be extremely small for 
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any realistic/natural symmetry breaking scale (for example, Planck, GUT or electroweak scales). 
Nonetheless, as a crude first approximation, we may assume that the formulae in Eq. (|4.1ip - 
(j4.12li hold true, even in the more general case, when the energy is emitted as a mix of massive 
and relativistic particles. 

There are a number of astrophysical observations and phenomena that can be interpreted in 
the framework of the cosmic string theory. Of these, gamma ray bursts (GRBs) [36] are one of 
the most extensively investigated. GRBs are powerful cosmic explosions that can be separated, 
based on their durations and characteristic wavelengths, in two classes m- short duration and 
hard spectrum, or long duration and soft spectrum, respectively. The possibility that GRBs are 
related to superconducting string emission via synchrotron radiation was investigated in detail 
in |38j-[42j. By assuming that the emitted energy is focused along a relatively narrow beam, 
the power released in a GRB event is of the order of 10^^ erg/s [36]. By assuming a similar 
beaming mechanism for the radiation emission from lumpy strings, it follows that a standard 
GRB type explosion can be generated by a single-lump string with A = ~ ~ 1 fm. On 

the other hand smaller cosmic lumpy strings can easily produce outputs of power of the order of 
10^3 _ ]^q 54 0 ];-g/g^ as observed, for example, in the (anomalous) cases of GRB 971214 and GRB 
990123 [21]. 


Lumpy strings may also be a powerful source of matter emission, with the emitted power 
almost of the same order of magnitude as the radiation power. It has been shown in |25j that 
a tangle of light, superconducting strings in the Milky Way could be the source of the ob¬ 
served 511 KeV emission from electron-positron annihilation in the Galactic bulge. In principle, 
powerful positron emission can also take place from lumpy cosmic strings. By assuming that 
some fraction r]^+ of the total energy is emitted in the form of positrons, we can write the re¬ 
lation (AE|„|)^_^ /At = {AN^+/At)jLmeC^ = Ve+ where AN^+ is the number of 

positrons, rUe is the electron/positron rest mass, and yp is the Lorentz factor of the positrons. 
Hence we can estimate the number of the positrons emitted per unit time by the lumpy cosmic 
strings as 


ANg+ 1 h |npa^ 

At 2^1 mee/. r/. 


(4.14) 


or 


AiVe+ 

At 


2.751 X 10^^ X ri^+ x 




I |2 2 

n a 


-1 


(4.15) 


The observed positron production rate in the galactic center is of the order of AN^^'^ ~ 1.2 x 10^^ 
positrons/s [25]. By assuming that 1% of the produced particles are in the form of positrons, 
ryg+ = 0 , 01 , the observed positron emission can be explained by the presence of lumpy strings 
at the galactic center, having rg = (again with |np 0{1)). Once produced by 

the decaying string lumps, the positrons will interact with the interstellar medium. The most 
common interaction processes are pair annihilation with an ambient electron, and positronium 
formation |25j . As a result of these interactions between positrons and interstellar matter, a 
narrow 511 keV line results, whose presence has been conhrmed by observations [25| . 


4.3 Cosmological applications 

However, cosmologically, perhaps the most important implication of the existence of lumpy 
strings would be their effect on the formation and evolution of string networks. In order to 
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investigate this, detailed models of the field dynamics during the symmetry breaking epoch must 
be developed in order to determine whether variations in local vortex radii may be correlated over 
regions A 2 Tc, as required for macroscopic string sections, filled with propagating lumps, to 
exist. In principle, this introduces a second “correlation length” into the network, in addition to 
the usual distance d « xct (x ^ I); which determines the average string length [21 [3]. Naively, we 
may expect this to signihcantly affect the growth of small scale structure for ^ d, smoothing 
out kinks over a time scale At ~ A^/c after the epoch of string formation, or to affect the 
macroscopic network dynamics for A^ 3> d. 

In the former case, the initial lumpy strings should evolve into smooth Nielsen-Olesen strings 
over a time. At, that is insufficient to affect the dynamics of the macroscopic network. In the 
latter, the effect of lumps propagating on the network is unknown. For small oscillations in 
the radial prohle, lumpy strings may behave just like uncharged current carrying strings, such 
as chiral strings [43|-[46] or neutral vortons m-m\ but, for large radial variations, they may 
resemble necklaces |12]-|18] with “beads” traveling at light speed. 

In general, lump propagation should damp traveling wave solutions, such as those presented 
in [30], when the lump width and amplitude are both smaller than those of the transverse mode, 
but it would be interesting to investigate the effect of resonance. It may be impossible to excite 
transverse oscillations, at a point z, whose amplitudes are smaller than the local vortex radius, 
or whose wavelengths are shorter than the lump width. Similarly, lump propagation may help 
to stabilize loops, of initial length p{ti) <C A 2 , over the characteristic time scale for lump decay, 
but momentum carrying loops of length p{ti) < A 2 could not form, and nor could loops with 
radii smaller than the maximum local string width. Given the possible complexity of all these 
interactions, it is difficult to estimate what overall effect the presence of lumps would have on 
the frequency of intersecting strings chopping of from the network, or the ability of the latter 
to reach a scaling solution [21E]. 

As an interesting toy model, but perhaps physically pathological case, one could investigate 
the possibility of local over-densities being sufficiently high to induce gravitational collapse, 
leading to the spontaneous production of black holes along the string. Accretion rates onto 
string networks |56]-[58]. gravitational radiation and lensing properties (see mm and references 
therein) and contributions to CMB anisotropies [59]-[ 66 ] should also be affected by lumps. 


5 Conclusions 


We have demonstrated the existence of a new solution of the abelian-Higgs field equations, which 
describes a generalization of the Nielsen-Olesen string. Physically, the new solution represents a 
non-cylindrically symmetric held conhguration, with variable scalar and vector core radii or, in 
other words, a lumpy cosmic string. In addition to the standard string parameters, such as the 
symmetry breaking energy scale p and topological winding number n, it is characterized by two 
arbitrary functions, obeying the standard relativistic wave equation, which represent arbitrary 
left or right movers. For a string with a static central axis, lying along the z-direction, these 
functions determine the local values of the t and ^-dependent core radii, so that the string carries 
momentum and the “lumps” propagate at the speed of light. 

Though symmetry breaking from the abelian-Higgs action remains a toy model for string 
formation, there is every reason to suspect that similar, generalized, solutions may be found for 
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more realistic string species, such as electroweak misiE], or GUT scale strings [7]. In effect, 
this introduces an additional phenomenological degree of freedom to the string network, which 
must be accounted for in any effective action, giving rise to position and time-dependent finite 
thickness corrections. In principle, there is no maximum value for the local string radius in the 
classical theory, though further analysis is required to determine whether a limit arises due to 
quantum effects. 

As such, lumpy strings embody several interesting features of existing models simultaneously. 
If the radial oscillations are highly localized, they resemble necklaces, with highly energetic beads 
traveling at light speed but, when individual lumps are highly dispersed along the string length, 
they may be seen as a form of neutral current, propagating in an (approximately) cylindrical, 
superconducting string. 

Since the lumpy string energy density is higher than that of the Nielsen-Olesen string, we 
expect the held conhguration to “relax” to Nielsen-Olesen form over some characteristic time 
scale, given approximately by the lump width divided by the propagation speed, i.e. the speed of 
light. With these simple assumptions, we were able to estimate the power emission from decaying 
lumps and found that they provide a viable mechanism for explaining anomalous gamma ray 
bursts |24] . in the case of radiation dominated emission, and may contribute to the observed 
positron excess emanating from the galactic centre [25], if the decay is dominated by matter 
emission. 
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